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ABSTRACT 

We consider sputtering of dust grains, believed to be formed in cooling supernovae ejecta, under the 
influence of reverse shocks. In the regime of self-similar evolution of reverse shocks, we can follow the 
evolution of ejecta density and temperature analytically as a function of time in different parts of the 
ejecta, and calculate the sputtering rate of graphite and silicate grains embedded in the ejecta as they 
encounter the reverse shock. Through analytic (ID) calculations, we find that a fraction of dust mass 
(1-20% for silicates and graphites) can be sputtered by reverse shocks, the fraction varying with the 
grain size distribution and the steepness of the density profile of the ejecta mass. It is expected that 
many more grains will get sputtered in the region between the forward and reverse shocks, so that our 
analytical results provide a lower limit to the destroyed fraction of dust mass. 

Subject headings: ISM : Supernova Remnants, ISM : Dust, Extinction, Galaxies : High-Redshift 



1. INTRODUCTION 

Understanding the chemical history of heavy elements 
( "metals" ) in the interstellar medium (ISM) of galaxies is 
an important, but difficult undertaking. Because refrac- 
tory elements are found to be depleted from the gas phase 
(locked into grains) , they provide a substantial reservoir of 
coolants in the solid-state phase. Interstellar grain surfaces 
are believed to furnish formation sites for molecular hydro- 
gen (H2), which is both a coolant and the starting point 
for a rich cloud chemistry. Grains also provide a critical 
transfer mechanism for reprocessing UV/0 starlight into 
far-infrared and sub-mm emission from interstellar molec- 
ular clouds and high-redshift galaxies. In the high-redshift 
universe, the first dust grains may influence the thermo- 
dynamics of the primordial gas and thereby control rates 
of star formation. 

Recent observations of dust grains at high redshift have 
been puzzling in the context of their formation. Several 
observations of damped Lya systems (Pettini et al. 1994; 
Ledoux, Bergeron & Petitjean 2002) have shown evidence 
of dust grains in them. Also, thermal emission (at 1.2 
mm) from dust in high redshift (z > 6) QSOs in the Sloan 
Digital Sky Survey (SDSS) have been reported (Bertoldi 
et al. 2003). The implied dust mass of 10^ Mq within 
a Gyr after the Big Bang appears difficult to explain with 
traditional models of grain formation in evolved low mass 
stars (Dunne et al. 2003), from which grains are thought to 
be transported to the interstellar medium (ISM) through 
stellar winds (Whittet 1992). 

Recent models have, therefore, focused on other sites of 
dust grain formation, such as ejecta of core-collapse su- 
pernova explosions that can occur on shorter time scales 
than that of the evolution of low mass stars (Dwek & Scalo 
1980). Core-collapse supernovae are thought to be capa- 
ble of synthesizing a significant amount of dust in the SN 
ejecta. Applying a theory of nucleation and grain growth 
developed by Kozasa & Hasegawa (1987), Kozasa et al. 



); 1991) and Todini & Ferrara (2001) have studied 
grain formation in expanding SN ejecta, the latter esti- 
mating the dust mass formed ~ 0.1-0.3 Mq. Nozawa et 
al. (2003) have extended it to the case of zero-metallicity 
SNe for population III stars (see also Schneider, Ferrara & 
Salvaterra 2004). 

Although the SN "dust factory model" would explain 
the large observed dust masses in high-redshift galaxies, 
observations of dust in nearby SNe appear inconsistent 
with this model, falling factors of 10-100 short of the re- 
quired amounts (~ O.2M0; see above). As an example, 
in a review of Type II SN 2003gd (NGC 628), Sugerman 
et al. (2006) state that "radiative transfer models show 
that up to 0.02 solar masses of dust has formed within 
the ejecta" . However, a recent study (Meikle et al. 2007) 
of the same remnant concludes that the mid-IR flux "is 
consistent with emission from 4 x 10^^ Mq of newly con- 
densed dust in the ejecta" . In a young, oxygen-rich SNR in 
the SMC, Stanimirovic et al. (2005) found only 10"^ Mq 
of hot dust {Td 120 K). The claimed sub-mm detection 
(Dunne et al. 2003) of ~ 3 Mq of cold dust toward the 
Gas A remnant was later shown to arise from interstellar 
dust in an adjacent molecular cloud (Krause et al. 2004). 

Observational tests of predictions made by models of 
SNe are difhcult, involving far-IR and sub-mm observa- 
tions of young SNRs. One must disentangle far-IR back- 
grounds, separate newly synthesized dust from circumstel- 
lar and interstellar dust, and understand the role of SNR 
reverse shocks, which can destroy newly synthesized dust 
in high-speed ejecta. A case in point is SN 1987A. Dust 
formation was inferred in ejecta of SN 1987A, but only 
10-3 Mq was detected (Dwek 2006; Dwek & Arendt 2007). 
Additional dust is almost certainly present, because of ob- 
served changes in the optical and bolometric fluxes and 
emission-line asymmetries in the red and blue wings (Mc- 
Cray 1993, 2007). Moreover, the reverse shock has not yet 
reached the interior of the supernova debris in SN 1987A. 
These inner ejecta are most prone to dust formation, since 
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they arc cold, dense, and metal-enriched. 

A challenging theoretical issue is whether high-velocity 
dust > 1000 km s^^) can survive the shocks, ther- 
mal sputtering, and grain-grain collisions, as the ejecta 
arc slowed down by the surrounding ISM. A goal of this 
paper is to determine how much dust survives the SN event 
and is incorporated into the surrounding ISM. The rela- 
tive speed of reverse shocks with respect to the expanding 
ejecta can be very large, and can raise the temperature 
of the ejecta to the extent of destroying the nascent dust 
grains in it. It is therefore important to study the effect 
of reverse shocks in detail. 

Recently, Nozawa et al. (2007) have considered this 
problem and came to the conclusion, with the help of hy- 
drodynamical simulation, that reverse shocks can destroy 
a fraction of dust mass of order 20-100% depending on 
the ambient gas density and explosion energy. Bianchi 
& Schneider (2007) also used numerical simulations and 
semi-analytical methods to conclude that only a small 
amount of dust mass survives; e.g., ~ 7% dust mass sur- 
vives in the case of a SN with a progenitor of mass 20 
Mq and ambient particle density of ~ 1 cm~^. Nozawa 
et al. (2007) have included the motion of dust relative 
to gas within the SN remnant and calculated the evolu- 
tion of dust grains over a period spanning the radiative 
phase (~ 10^-10^ yr), whereas Bianchi & Schneider (2007) 
neglected the effect of motion of grains on dust destruc- 
tion and calculated dust evolution up to the non-radiative 
phase (< 10^ yr). 

We study in this paper a simple and analytically 
tractable problem, namely the effect of reverse shocks for a 
power-law mass distribution of the ejecta, in the regime of 
self-similar evolution of both forward and reverse shocks. 
Although the simple formalism does not allow us to study 
several important processes — including motion of dust 
grains and its effect on evolution of dust grains — our ap- 
proach allows one to segregate the processes of destruction 
within the reverse shock itself from those occuring between 
forward and reverse shocks (which can be dealt with in nu- 
merical simulations involving grain motion). 

2. DUST IN SUPERNOVA REMNANTS 

We first summarize the results of the theoretical pre- 
dictions of dust grain formation in SNRs. Todini & Fer- 

rara (2001) calculated the formation of different types of 
dust grains — amorphous carbon (AC), silicates (enstatites 
: MgSiOz, forsterites : Mg^SiO^) and iron-bearing mag- 
netites {Fe^Oi) and others, in ejecta with different metal- 
licities. They found that AC grains are generally of larger 
size, in the range of ^ 0.01 0.1 /i m, whereas silicates are 
produced with sizes ^ 10^'^ m. Nozawa et al. (2003) 
reached similar conclusions for Pop III SN events, with 
sizes of newly synthesized grains spanning a range of three 
orders of magnitude, with maximum grain size being less 
than ^ 1 /im. 

The dust grain parameters in the ISM are inferred to be 
somewhat different from these estimates. The extinction 
observations are satisfactorily reproduced by assuming 
mostly of two components — graphite and silicate — with a 
common power-law size distribution, dn/da oc a~^'^, trun- 
cated at a minimum size ai~5xl0~^ jim and a maximum 
size am ^ 0.25 /xm (Draine 2003, and references therein). 



There also appears to be a population of very small grains 
with PAH composition. Nozawa et al. (2003) found for 
grains synthesized in Pop III SN ejecta a size distribu- 
tion with index —2.5 for smaller grains and with —3.5 for 
larger grains (with the grain size at crossover point de- 
pending on supernovae models). We use these two power 
laws to bracket the possible size distributions of grains syn- 
thesized in SN ejecta, and assume that grains are formed 
in the range 10~''~3 x 10~^ cm (see Figure 11 of Nozawa 
et al. (2003). Wc also consider the; c;ase of an upper limit 
of 10~^ cm and compare our results for both upper limits. 

It is instructive to estimate the destruction time scales 
before we calculate the effect of reverse shock in de- 
tail. The sputtering time scale for dust (graphite, sil- 
icate and iron) grains of size a is approximately ~ 
^^^^ (r^) ( I yr, for ambient temperatures larger 

than 10^ K (Draine & Salpeter 1979), where Ue is elec- 
tron density. For grains with size ~ 0.1-0.3 /xm (which 
incidentally contain most of the dust mass for a a~^'^ dis- 
tribution), and for Ug ^ 1-10 cm~^, the sputtering time 
scale is ~ 2 X 10^~^ yr. 

3. REVERSE SHOCKS 

We now briefly discuss the self-similar evolution of for- 
ward and reverse shocks in SN (see Truelove & McKee 
1999, hereafter referred to as TM99, for further details). 

3.1. The ejecta 

The parameters used to describe a supernova remnant 

which affect its evolution arc the energy of the explosion 
(Eej), the mass contained in the expanding ejecta (Mej), 
the maximum velocity of the material within the ejecta 
{vej), and the density of the ISM (po)- Denoting the time 
of the explosion as t = 0, we construct our initial con- 
ditions at a later time, still early in the history of the 
remnant, assuming the ejecta to have already expanded 
to a radius Rej without any deceleration due to the am- 
bient medium. If we consider SN to be a point explosion, 
then the radial velocity profile will be linear as ejecta with 
velocity v will travel out to a distance vt : 

This profile becomes flatter as time progresses and the 
faster-moving ejecta arc flung out to greater distances. 
The corresponding radial density distribution of the ejecta 
matter is given by: 

p{r) = < '^"^^ ' ' K^J^^^3' ' (2) 

[ po, T> Rej 

where f{v/v,,j) is the time-independent structure fimction 
of the density profile and the term arises from the free 
expansion of the ejecta. We define a dimensionless param- 
eter to label shells of the ejecta starting from the centre 
outwards as 

w= —,0<w <1 (3) 

Following TM99 we consider the situation when / is given 
by a power law in w: 

f{w) = I _„ « ^ (4) 

^ ^ \ fnW Wc<W<l ^ ^ 
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where we have separated the ejecta into a uniform density 
core region and a power-law envelope region. The conti- 
nuity of / at Wc and the normalization of p translate to 
expressions for parameters /o and /„ in terms of the free 
parameters Wc and n: 

/O = fnW-^ 



fn — 



3 

4-77 



n/3 



1 - in/3)wi 



(5) 



3.2. Forward 



ana reverse 



shocks 



As the freely expanding ejecta come into contact with 
the ambient medium, with the contact discontinuity mov- 
ing at Vej, a speed much larger than the sound speed in 
the ambient medium, a forward shock is set up in the 
surrounding gas which propagates into the medium. The 
ejecta at the contact surface decelerates suddenly as it hits 
the ambient gas, and consequently, a reverse shock is set 
up that propagates inwards through the ejecta. While the 
mass of this matter swept by the forward shock is less than 
the total mass of the ejecta driving the shock, the remnant 
is said to be in the Ejecta Dominated or ED stage. It is 
in this stage, before the reverse shock has attained a sig- 
nificant velocity, that the ejecta can be approximated as 
a spherical piston freely expanding into the ambient gas. 
When the mass swept by the forward shock becomes ap- 
proximately equal to the ejecta mass, the corresponding 
late-time limit is known as the Sedov-Taylor or ST stage 
of the evolution of the remnant. 

Wc follow TM99 to follow the evolution of the reverse 
shock analytically. Let Rj and i?,. denote the radii of 
the forward and reverse shocks, and Vf = dRf/dt and 
Vr = dRr/dt be the velocities of the forward and reverse 
shock waves respectively, all in the rest frame of the un- 
shockcd ambient medium. We have the velocity of the 
unshocked ejecta as seen by the reverse shock as it propa- 
gates into the ejecta: 

R (t) 

Vr = v{Rr, t) — Vr{t) = Vr{t) . (6) 

Let (f) be the ratio of the pressures behind the reverse and 
forward shock waves and I be the ratio of the radii of the 
two shocks at a given time: 



pav}{t) 



Rr{t) 



Wf{t) 



Rf{t)/t 



It is useful to express various parameters in dimension- 
less forms in the units of characteristic values of these vari- 
ables, so that, Rt, = R/Rch, t* = t/tch, and = v/vch- 
For example, the characteristic scales of length, time and 
velocity are given by. 



Rr 



l/3„-l/3 



M^ypo 



tr.,=E-'/'M!/\ 



5.25 n, 

-1/3 



-1/3 







pc. 



1.66 X Wn, 



-1/3 







Vch = Rch/tch « 3.16 X 10^ km s-i . (8) 

(We have assumed po = lAnorUp for He/H= 0.1 by num- 
ber.) We have estimated the numerical values for the 
case E = 10^^ erg, Mej = 10^^ g and no = 1 cm^^. 
We wiU therefore write R* = R/Rch, t* = t/tch, and 
V* = v/vch, and the radius of the reverse shock, in partic- 
ular, as R* = Rr/Rch- 



3.2.1. ED (Expansion dominated) stage 

The values of (f) and I approach constant values as t ^ 
if the solutions remain self-similar in the same limit. We 
write (f){t) ~ 0(0) = and l{t) ~ /(O) s and also 
w/(0) = I. The solution to reverse shock evolution ex- 
ists in the form of two branches, which we label the en- 
velope {wc < Wf/lED < 1) and core (0 < Wf/lED < Wc) 
branches. Denoting a = ... — 2-, the evolution of the 

reverse shock is given by, 

t*{R;) = {^yR;x 

-2/(3-n) 



1 



Wc 



*3/2 



2/(P) 



where the first part corresponds to the case < Wf /I < 1 
and the second part, to < Wf/l < Wc- The corresponding 
velocity of reverse shock in these two cases are, 



v*{r;) 



1/2 



5/3 



(10) 



The relative velocity of the reverse shock in the frame of 

reference of the expanding ejecta, for the case when the 
reverse shock is in the envelope region, is given by. 



I^ir;'/' 



( 




2/(3-n) 


\ 






(11) 



describes the reverse 



Note that, since Wf/l 

shock, we will define it as and use it to label shells 
in the ejecta. 

This solution is valid (in the ED stage) only as long as 
the reverse shock remains in the envelope Wcore < < 1- 
In the case of n < 3 ejecta, there is no need for a core and 
wc can take the limit Wcore ^ 0. For n > 3 however, we 
must have a "flat profile" core to keep the integrated mass 
finite. For the steep index scenario, therefore, the func- 
tional form of the solution will change when the reverse 
shock reaches the core and we will return to this point 
later. 

3.2.2. ST (Sedov-Taylor) stage 

We now describe the Sedov-Taylor stage, which is the 
late-time evolution limit of the non-radiative phase of a 
supernova remnant. Following TM99, we assume that the 
shocks trace a constant acceleration path after transition 
to the ST stage at tsT, which (approximately) begins when 
the mass of the ambient medium shocked by the forward 
shock equals the mass ejected by the supernova. Thus the 
velocity of the oncoming ejecta in the frame of the reverse 
shock is 

K{t*)=K,ST + K,ST{t* -^*St), (12) 
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where hr^sT is the constant acceleration and Vr^ST is the 
velocity of the ejecta w.r.t the reverse shock at tsx- Note 
that Vr = —td{Rr/t)/dt. Using this and integrating the 
resulting expression from tsx to t gives 
'R 



Rl(t*) = t* 



ST 



a. 



■r,ST''ST 



ST 




For the n < 3 case, we therefore have an implicit solu- 
tion for the reverse shock position in the ED stage t < tsT 
and an explicit solution in the form of a uniform acceler- 
ation description in the ST stage. For the n > 5 case, we 
again have an implicit solution in the ED stage at least 
while t < tcore- In fact, an explicit solution for R*{t*) can 
be determined from the implicit solution in the ED stage 
for n > 5 remnants. We have (see equation (74) of TM99), 

1/2 1 -2/(3-n) 



r(i?:)^(-) r; 



In this expression, we can take the limit Wc 



(14) 
by 



allowing Ve 



oo, the ejecta energy E remaining finite. 



In this limit fn ^ and when the second term in the 
brackets becomes much larger than unity, we have 

27 1 flO /n-5 



R;{t*) 



4Trn{n-3)P(l) 
This gives us the velocity of the reverse shock as. 



(15) 



27 



1 



47r n(n - 3)^2,/, 



10 

y 



n — 5 
n-3 



^1 1/ 



and correspondingly, 

v^{t ) = 



(16) 



t 



27 



1 



47r n(n - 3)P(t) 



10 

y 



-, l/n 



t*-li7) 



This is the result arrived at by Chevalier (1982) and Nady- 
ozhin (1985) - henceforth referred to as 'the CN solution'. 
The transition of a supernova remnant with n > 5 to this 
stage has been shown to occur extremely fast (TM99). 
Thus for all practical purposes, we may assume that the 
remnant begins with the CN stage. 

The reverse shock then moves into the core at t = tcore- 
TM99 assumed that the reverse shock motion after this 
time is also described by a Sedov- Taylor trajectory - which 
leads to better match between the two regimes in terms of 
continuity of the shock radius and velocity. In the case of 
n > 5, therefore, if we replace tsT by tcore in the ST stage 
equations above, the same analysis follows through. The 



numerical values of the parameters t, 



core 5 '^TjCore 



and Vr 



are, of course, different from the parameters for the n < 3 
case: viz tsT and ar.ST- We have, 

_ D + 



R;{t*)^t* 



t* 



.it* 



t 



^r, core core 



)ln 



Kin 



(t* - t* ] 



(18) 



The values of the parameters used in our calculations have 
been taken from TM99, which have been checked with re- 
sults from simulations. 

3.2.3. Summary of ED and ST solutions 

lel thus have the final solutions for the reverse shock 
as follows. In the case of n < 3, equations ([U first part) 
and (fTTjl describe the solutions in the ED stage, and, there- 
after, equations and are appropriate for solutions 
in the ST stage. For the case of n > 5, solutions are given 
by equations and ([TT]) in the t < tcore stage. Af- 
terwards, in the t > tcore stage, equations (fTS]) describe 
the solutions. The analytical solutions remain valid until 
t^ ~ 2.2 (TM99), and we calculate the evolution of shocks 
until this point (except for n = 4; see below). 

For a given set of values of E and Mej , we can therefore 
solve for reverse shock position and relative velocity, as- 
suming a density distribution described by index n. The 
parameter a in different cases of n, is given by. 



E 



{l/2)Mc,v% 



(5-n)( 



-(5- 

'Wcore 



n/5 



-(3- 
UJcore 



n/3 



For n < 3, we calculate this in the limit of Wcore = 0, 
and for n > 5 (including the case of n = 4) we assume 
(following TM99) a value of Wcore = 0.1. Values of other 
parameters (such as are taken from appropriate tables 
of TM99 for the relevant values of n. 

We also consider the special case of n = 4, for which 
the mass integral diverges without a core. We follow 
-^%tee prescriptions of TM99 (their section 8) for this case 
(and, again, assume Wcore = 0.1), and track the reverse 
shock with only ED solution without any transition to the 
ST case. According to TM99, based on their compari- 
son of analytical and numerical results, this prescription 
gives tolerably good results (accurate to within 24%) until 
t* = 1.2. 

3.3. Arrival time of reverse shock 

We have introduced a label to a given shell in the super- 
nova ejecta hy w = Wf /I = ^''J*^^* , the ratio of the velocity 
of the ejecta in the shell to the maximum ejecta velocity. 
The ejecta in a shell specified by a given w first undergoes 
free expansion (in the ED stage) and self-similar expansion 
(in the ST stage) until hit by the incoming reverse shock 
wave. This occurs at a time, which we call thu, which is 
a function of w and which is found by equating w to the 
Vr/vej found from the position of the reverse shock. That 
is, 

Vr Rr ^ jy jy 

w = — — — — ^ Rr ~ wR 

Vej ^ej 

Or, in terms of dimensionless variables. 



WVedhit 



(20) 



Rl = w ■ 



Vejthit 

R. 



■ch 



Veit, 



ej '^ch 



R. 



■ch 



hit ' 



(21) 



This expression for R* in terms of t* is equated to the 
explicit solution for the reverse shock, when available, as 
a function of time, and the resultant linear equation in 
tj^jj and w is solved using the bisection method for each 
given w. This is applicable to the steep (n > 5) ejecta case 
(both the CN and ST stage) as well as to the ST stage for 
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shallow {n < 3) ejecta. For the ED stage of shallow ejecta, 
we use t^jj as a function of R* in equation ([9]) and then 
employ the bisection method exactly as before to solve for 
^hif Either way, we can determine i^^j given any shell 
labeled by a w in the SN ejecta. 

3.4. Density and temperature evolution in ejecta shells 

We assume that gas in a given shell of the ejecta evolves 
adiabatically before thu- The temperature is then raised 
to the postshock temperature and the density acquires a 
jump, corresponding to the value of relative shock velocity 
Vr- After thiti the gas again evolves adiabatically, that is, 
its density decreases as t~^ (equation [J). Its temperature 
decreases as T cc p''~^ cx for 7 = 5/3 appropriate for 
monoatomic gas. For a shell with a certain value of w, 
therefore, we have the following evolution of density. 



pit,w) 



and temperature 
T{t,w)-^ ^'^^^ 



T, 



t ^-2 

r shock ( . 



t < thit 
thit < t < '2.-2tr_h 

t < thit 

thit < t < 2.2tch 



(22) 



(23) 



where Trshock — j^fJ-'TipVr'^ /ItB, the post-shock temper- 
ature due to the reverse shock with relative velocity Vr- 
For a completely ionized gas the mean molecular weight 
^ = 0.6. We assume a value of T; = 5400 K at = 70 
days, motivated by photometric observations of SN1987A 
(Catchpole et al. 1987) where it was estimated that the 
photospheric temperature at ~ 70 d after the explosion 
was ~ 5400 K. We have also checked that our results do 
not depend strongly on these assumptions. 

3.5. Sputtering of dust 

Given the above mentioned evolution of density and 
temperature in shells of ejecta labeled by different val- 
ues of w, we can compute the steady decrease in grain 
radius, for different types of dust grains. We use the poly- 
nomial fit given by Tielens et al (1994) for graphites and 
silicates as a function of gas temperature. They expressed 
the rate, (I/uh)^ in powers of logj^Q T, with coefficients 
for different grain composition tabulated in their Table 4 
(see their equation 4.21). These rates are consistent with 
recent results of Nozawa et al. (2006). 

We follow a shell of ejecta material labeled by a value 
< w < 1 from ~ 10 years after the explosion (by which 
time dust grains are believed to have formed in the cool- 
ing ejecta, since the formation time scale is of order a few 
years (see, e.g., Todini & Ferrara 2001)), to when the re- 
verse shock hits the shell at t^it and beyond, until a time 
t* ^ 2.2, which is the limit of the validity of the analytical 
solutions. 

4. RESULTS 

We have calculated the effect of reverse shocks on sput- 
tering of graphite and silicate grains for typical value of 
explosion energy, stellar ejecta mass and ambient density. 
Below we show the results for E = 10^^ erg, Mej = lO"^** 
g, po = lO"^'* g cm"'\ 

First, we show the evolution of grain radius (for 
graphites and silicates) for two dust grains embedded in 



two particular shells, for example, in ejecta characterized 
by n = and E = 10^^ erg. We show in Figure 1 the cases 
for the shells marked w = 0.3, and 0.6. The reverse shock 
passes through these two shells, first hitting the w = 0.6 

shell at ~ OMch ~ 0.8 x lO'^n^^^^ yr, and then the w = 0.3 
shell at l-2tch ~ 2 x lO'^n^ yr. The bottom panel of 
Figure 1 shows the rise in temperature for these two shells 
followed by adiabatic cooling. Initially the temperature 
drops to very low temperatures due to strong adiabatic 
cooling, although, in reality, gas cannot be cooled below 
the temperature of the cosmic microwave background, but 
we have not used such a lower bound in our calculations. 
At any rate, the post-shock temperature does not depend 
on the low temperature of the gas upstream, being deter- 
mined only by the strength of the reverse shock, and so 
the results of our calculation remains realistic in spite of 
the strong cooling before the reverse shock hits the gas in 
a given shell. 

Figure 2 plots the ratio of final to initial grain sizes as a 
function of the initial grain sizes for silicates and graphites 
forming (and being sputtered) in shells w = 0.2, 0.4, 0.6 for 
the particular case of n = 0. Since the sputtering yield of 
silicate grains is larger than that of graphites for a given 
set of parameters, silicates are eroded more rapidly than 
graphite grains for sputtering in a given ejecta shell. Also, 
the relatively outer shell oi w = 0.6 experiences less sput- 
tering than the inner shells oi w = 0.2 and 0.4. This is 
expected since the reverse shock (relative) speed picks up 
as it plows through the ejecta, increasing with decreasing 
values of w. To illustrate this point, we plot the trend of 
the relative velocity of the reverse shock as a function of 
w for the shells hit by the shock in Figure 3, for n = 
(solid), n — 2 (dotted), and n — 6 (dashed). The relative 
velocity is a rising function with decreasing w, for n = 
and for n > 5, and it is a peaked function for n = 2 (and 
also for n = 4, which is not shown here, but as shown by 
TM99). 

However, we notice in Figure 2 that grains (both 
graphites and silicates) in shell w = 0.2 (solid curve) are 
sputtered to a somewhat lesser extent than grains in the 
w — 0.4 shell, although the relative speed of the reverse 
shock should be (according to Figure 3, solid curve for 
n — Q) higher when it hits the w — 0.2 shell than at the 
w — 0.4 shell. Although the post-shock temperature at 
w = 0.2 is higher than at w = 0.4, the shells in the inner 
region encounter the reverse shock later than shells in the 
outer region of ejecta, and by the time they encounter it, 
the density in the inner shells would have decreased rapidly 
(see equation . The rarefaction in the inner shells at 
the time of encounter with the reverse shocks lessens the 
effect of dust sputtering in these shells. 

We illustrate this point by considering grains which are 
sputtered below a certain size, say, 10""^ cm (10~^ /^m), as 
a function of the shell in which they reside (and get sput- 
tered). We determine the maximum size of grains which 
are sputtered below a size ~ 10^^ cm, as a function of 
w, for different shells, and present the results in Figure 4. 
Curves for different values of n show the maximum size of 
grains sputtered below 0.001 /^m as a function of the pa- 
rameter w (left panels), as well as against the fraction of 
ejecta mass that is contained inside of the concerned shell, 
M{< w)/Mej in the right panels. Results for graphites 
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and silicate grains are shown in the upper and lower pan- 
els respectively. 

Figure 4 clearly shows that sputtering becomes intense 
in intermediate shells, and decreases (as a result of rar- 
efaction) with decreasing value of w for small values of w. 
The sizes of grains which are sputtered below 10~^yum lie 
in the range of < 10~^ /im -10~^'^ fim, with silicates more 
strongly sputtered than graphites. In all cases, the curves 
of the maximum size of sputtered grains increases with de- 
creasing values oi w and M{< w)/Mej, finally plunging to 
very low values for w <C 1, indicating negligible effect of 
the reverse shock deep inside the ejecta. For example, the 
(dotted) curve for n = 2 shows a peak at w ^ 0.17 below 
which the maximum size of sputtered grains decreases. 

Also, although it is not explicitly shown in Figure 4, the 
corresponding results for n > 5 show that reverse shocks 
do not penetrate most of the ejecta mass in the time period 
of the validity range of analytical solutions ^ 2.2. Ana- 
lytical means allow us to probe deep into the ejecta mass 
distribution only for n < 3, and full hydrodynamical mod- 
els are needed to address the cases of steeper profiles. For 
steep density distributions, the limiting time of t = 2.2tch 
is reached by the time the reverse shock only skims the sur- 
face of the ejecta, as far as mass fraction is concerned. The 
characteristic time in the case considered here (E = 10^^ 
erg, Mgj = 10^^ g, ambient density po = 10~ g cm~^) 
is tch ~ 1.7 X 10^ yr. Although, by the limiting time of 
2.2tch ^ 3.7 X 10"^ yr for the validity of self-similar regime, 
reverse shocks can hit shells with small w in these steep 
cases, the mass contained inside of these shells remains 
very large. In other words, in these cases wc cannot probe 
what happens to most of the ejecta mass, or the dust con- 
tained in them with analytical means, within the limits of 
the self-similar regime of shock evolution. 

We also note that a substantial amount of grain sput- 
tering takes place in the hot gas between the forward and 
reverse shock that analytical methods cannot capture. De- 
tailed hydrodynamical calculations are needed to address 
these issues, and our result can only provide a lower limit 
to the total mass fraction of dust that is destroyed. 

Finally, wc calculate the total dust mass sputtered away 
as a function of initial grain sizes in the following manner. 
We first calculate for a given shell w, the final grain size 
af{ai,w) corresponding to the initial grain size aj. Then 
the fraction 1 — {af/ui)^ is the total dust mass sputtered 
for grains formed in this shell. We then calculate the to- 
tal dust mass sputtered for all shells (for a given initial 
grain size Oj), weighted by the shell mass in which they 
are situated, as given by, 

where M(< w) is the ejecta mass contained within the 
shell w. 

We plot these weighted fractions fd{ai) for graphite and 
silicate grains as functions of for n 0,2,4 in Figure 
5. Thick curves denote the results for silicates and thin 
ones are for graphite grains. As expected, the sputtered 
fraction decreases with increasing a^, and the fraction for 
silicates is in general larger than that for graphites. 

The total dust mass fraction that is destroyed can then 
be evaluated from these results by convolving fd{ai) with 
an initial grain size distribution. As mentioned earlier, 



Nozawa et al. (2003) found that the size distribution of 
synthesized grains is bounded by two power laws, with in- 
dices —3.5 for larger grains and —2.5 for smaller grains. We 
therefore estimate the total dust mass fraction destroyed, 
assuming a size distribution n{a)da oc a^^da, with p = 2.5 
and 3.5. The destroyed fraction of dust mass is given by, 

^ ^ J n{a)a^ fd{a)da 

J n{a)a^da ' 

where the grain volume scales as a^. The results for 
graphites and silicates are presented in Table 1, for n = 
0,2,4 {E = 10^1 erg, M^j = lO^"* g), leaving the cases of 
n > 5 for which the reverse shocks do not reach the inte- 
riors of the ejecta within the validity range of self-similar 
solutions. We assume the lower and upper limit for the 
integral to be 10~^ cm and 3 x 10~^ cm respectively, as 
explained earlier in §2. We also calculate the destroyed 
fractions for a truncated grain size distribution, with an 
upper limit of 10~^ cm, and the results are shown inside 
brackets in Table 1 for comparison. 

We have further studied the variation of the destroyed 
dust mass fraction with explosion energy and ejected mass. 
Tables 2 and 3 show the fraction of dust mass supttered 
for two sets of parameters: Table 2 for M^j = 10'^'* g and 
E = 10^2 erg, and Table 3 for M^j = 2 x lO^" g and 
E = 10^^ erg. Figure 6 plots the varition of (logarithm) of 
the destroyed dust mass fraction, /, for silicates, for dif- 
ferent values of explosion energy and ejected mass, and for 
different mass profiles {n = 0, 2, 4). In the upper panel, we 
plot the variation of / with M^j (in logarithm), keeping 
the explosion energy constant dX E — 10^^ erg, and in the 
lower palnel, we plot (logarithms of) / with explosion en- 
ergy E, keeping the ejected mass fixed at Mej l034 g. The 
destroyed fraction generally increases with ejected mass, 
but its variation with explosion energy is not monotonic. 
The reason is that the process of dust destruction depends 
mainly on two factors: the speed of the reverse shock 
(which depends on the characteristic speed Vch) and the 
characteristic time, tch- Now, Vch oc M~^^'^ E^l"^ , whereas 

tch, oc mI'.'^E-^I'^. It would appear that a change in 
ejected mass manifests itself in increasing tch-, and thereby 
increasing the destroyed fraction, whereas, a change in E 
affects tch and Vch in a complicated manner, which results 
in the variation shown in the bottom panel of Figure 6. 

5. DISCUSSION 

We therefore find that a mass fraction < 20 % of silicates 
and graphite grains created in cooling ejecta can be sput- 
tered away by the reverse shock if the density distribution 
of the ejecta has a shalow profile (n < 5). Nozawa et al 
(2007) recently studied sputtering of dust in reverse shock 
and found a destroyed dust mass fraction of 0.2-1; their 
result for an ambient hydrogen number density of 1 cm~^ 
and a progenitor mass of 13 M© is ~ 0.7 (their Table 1). In 
fact, their result for the 'unmixed grains' model (which has 
to do with spatial mixing of the core heavy-element distri- 
bution) in which the maximum size is 0.01 /zm is ~ 0.5. It 
shoiild be noted, however, that a considerable amount of 
sputtering in their calculation is caused by the hot plasma 
between the forward and the reverse shocks (see §3.3 in 
Nozawa et al (2007)), which we have not considered. In 
our case, for analytical simplicity, we have assumed the 
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post-reverse shocked gas to be coohng adiabatically. Our 
results are, therefore, expected to give lower bounds on 
the actual dust mass fraction destroyed in SN remnants. 

There is another reason why the destroyed fraction of 
dust mass calculated here is lower than that of Nozawa et 
al (2007) and Bianchi & Schneider (2007). Strictly speak- 
ing, the destruction rate of Tielens et al (1994) used here is 
valid for a gas with solar abundance. The gas in the dust 
formation region in the supernova ejecta is expected to be 
metal rich, and the destruction of dust grains in this gas 
can be more efficient in this case, since sputtering yields 
by metal ions are much higher in general than by hydrogen 
and helium ions. 

It should be noted that our results for the destroyed 
fraction is independent of the ambient density, since the 
gas density enters into our calculation only to change tch 
(notice in equations [S] that Vch is independent of ambient 
density, and so the relative speed of reverse shock, and 
in turn the post-shock temperature, are also independent 
of ambient density). In reality, however, a lot of grain 
sputtering would take place between the forward and re- 
verse shock, where the gas density will crucially depend 
on the ambient density. This aspect of grain sputtering is 
neglected in our analytical calculation, but is captured in 
the simulation results of Bianchi & Schneider (2007) and 
Nozawa et al. (2007). 

Another difference between our results and those of 
Bianchi & Schneider (2007) arises from the assumption of 
initial grain size distribution. Bianchi & Schneider (2007) 
used grains of small size in their calculation, with silicate 
grains smaller than < 0.01 /im and graphites smaller than 
< 0.05 /im. Also, their size distribution is more confined to 
a narrow range than the power law distribution assumed 
here. Both these factors would enhance the destroyed 
fraction of dust mass in reverse shocks. For example, if 
the sizes of Si02 grains have a delta-function distribution 
around ^ 25 A (see Figure 1 of Bianchi & Schneider 
2007), then the curves in Figure 2 of the present paper 
(say, for w = 0.2) yield a mass destruction fraction of 
- 22 % (with Qf/a, - 0.92). For MgSiOs grains of size 
~ 45 A, we similarly estimate a destroyed mass fraction 
of ^ 14 %. From Figure 2, we estimate that one requires 
aj ^ 10 A in order to get ~ 50% destruction of silicate 
grains. 

A note on the observed density distribution of super- 
novae ejecta is in order here. Chevalier & Fransson (1994) 
studied a model of the ejecta from a Type II supernova 
with a shallow (n < 3) core surrounded by a steep {n > 5) 
envelope. In the case of Type la supernovae, a steep den- 
sity distribution in the ejecta mass has also been discussed 
in the literature (Colgate & McKee (1969); Dwarakadas 
& Chevalier (1997)). There has been success, however, 
in modelling data of Type la supernova with uniformly 
distributed (n ~ 0) ejecta (Hamilton & Sarazin (1984); 
Hamilton, Sarazin & Szymkowiak (1986)). 

It is worth considering what observations might improve 
our understanding of the net dust production by SNe. 
Far-infrared and sub-millimeter observations of reverse- 
shocked dust in young SNe by Spitzer, Herschel, and other 
facilities will help, but our theoretical studies suggest that 
unshocked (cold) dust may harbor a considerable mass of 
undetected dust. To constrain this dust and distinguish it 



from surrounding interstellar clouds will probably require 
sub-mm observations with small beam sizes. This brings 
in complications of the circumstellar environment of young 
SNe. Furthermore, the long times (thousands of years) re- 
quired for reverse shocks to reach the core ejecta suggest 
that late-time observations of supernova remnants would 
be useful. 

6. SUMMARY 

We have studied the effect of reverse shocks analytically, 
in the regime of self-similar evolution, on the sputtering of 
dust grains that are believed to form in cooling ejecta. 
For representative cases we found that fractions of dust 
mass that is destroyed are of order of order 1-20 % of 
silicates and graphites, for SN of explosion energy ^ 10^^ 
erg, ejecta mass ^ lO'^^ g, and an ambient density ~ 10"^'* 
g cm~'^. Our analytical formalism provides only a lower 
bound on the dust mass fraction that is destroyed by the 
reverse shock, since it ignores further sputtering of grains 
in hot plasma between the forward and reverse shocks. 
Our results are, therefore, consistent with the recent es- 
timates from numerical simulations which include these 
additional effects. Furthermore, Our study provides a for- 
malism to generalize these results to cases with different 
values of SN parameters. 

We thank Drs. R. Chevalier, E. Dwek, T. Nozawa, 
A. Ray and Y. Shchekinov for helpful discussions. This 
work was supported at the University of Colorado by 
NASA Theory grant NNX07-AG77G and NSF theory 
grant AST07-07474. 
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Table 1 : Fraction of dust mass sputtered for sizes ranging from 10 /um to Umax = 0.3 /xm (or Umax = 0.01 /um) 



grain (dn /(la) 


II = 


/( = 2 


= 4 


Graphite (oc a'"^-^) 


0.004(0.02) 


0.01 (0.05) 


0.007 (0.04) 


(oc a-^-^) 


6 X lO"'' (0.01) 


0.001 (0.03) 


0.001 (0.03) 


SiUcate (oc a""^ *^) 


0.01 (0.07) 


0.03 (0.13) 


0.02 (0.10) 


(oc a-^-^) 


0.002 (0.05) 


0.005 (0.10) 


0.004 (0.07) 



Table 2 : Fraction of dust mass sputtered for Mgj = 10^'' g, Egj = 10^^ erg. 



grain {jlii /da) 


;/ = 


II = 2 


;) = 4 


Graphite (oc a^^ '^) 


0.001 (0.008) 


0.003(0.01) 


0.005(0.03) 


(oc a-"-^) 


2 X 10-"^ (0.005) 


4 X lO"'' (0.01) 


7 X 10-'' (0.02) 


Silicate (oc a"^ *") 


0.006 (0.03) 


0.01 (0.06) 


0.05 (0.16) 


(oc a-^-'') 


0.001 (0.02) 


0.002 (0.04) 


0.018 (0.14) 



Table 3 : Fraction of dust mass sputtered for Mej = 2 x 10^** g, E^j = 10^^ erg. 



griiiu ((/;) /da) 


;( = 


;/ = 2 


;) = 4 


Graphite (oc a^'^ '') 


0.006 (0.035) 


0.015(0.08) 


0.01 (0.06) 


(oc a-'-^) 


10-^ (0.02) 


0.003 (0.06) 


0.002 (0.04) 


Silicate (oc a"^ *') 


0.02 (0.1) 


0.04 (0.18) 


0.03 (0.13) 


(oc a-^-'^) 


0.003 (0.07) 


0.008 (0.14) 


0.005 (0.10) 
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Fig. 1. — The dust grain size as a function of time (where = t/t^f^) is plotted for ii = 0, E = lO''^ erg, M^j = 10'^"^ g in the upper panel, 
for two shells, with w = 0.6 and w = 0.3. The shell w = 0.3 is inside oi w = 0.6 shell, and is hit by the reverse shock at a later time. The 
initial grain size (solid line for graphites, and dotted, for silicates) is 1.3 x 10~^ /im. The bottom panel plots gas temperature as a function 

of t* showing the rise at the instant of the reverse shock meeting the two particulaj: shells, {tch 
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Fig. 2. — Ratio of final to initial grain radius is plotted against initial grain radius (in /im) for silicates (thick lines) and graphites (thin 
lines), for grains in shells w = 0.2 (solid), w = 0.4 (dotted), w = 0.6 (dashed) for the case of n = 0. 
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Fig. 4. — The values of maximum size, \ogamax (in below which grains are sputtered below 10^^ cm are shown against the parameter 

w in the left panels, and against the fraction of ejecta mass inside the corresponding shell, M(< w)/M^j , in the right panels, (w = R{t) / (v^jt), 
or equivalently, it is the ratio of the velocity of the ejecta in the shell to the maximum ejecta velocity.) The case of graphites are shown in 
upper panels and the case of silicates, in the lower panels. Dust is introduced at 10 yr, where the tracks begin. Solid, dotted, and dashed lines 
refer to the cases n = 0, 2, 4 respectively. The energy of explosion is assumed to be E^j = 10®^ erg and ejecta mass is assumed to be lO^"* g. 
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Fig. 5. — Fraction of dust mass sputtered as a function of initial grain size, weighted by the mass of the shell in which they are formed, 
for n = (solid), 2 (dotted) and 4 (dashed), for silicates (thick curves) and graphites (thin curves). 



Dust sputtering by reverse shock in SN 



15 




34 34.5 
Log 10 Ki 




34 34.5 
Log 10 Ki 




34 34.5 
Log 10 



-1 



-1.5 



-2.5 



-3 



T — I — I — r 




: n=0 



J I I L 



51 52 
Log 10 E 




51 52 
Log 10 E 



-1 



T — I — I — r 




51 52 
Logio E 



Fig. 6. — Fraction of dust mass sputtered is plotted in the upper panel as a function of mass of ejecta, keeping the explosion energy 
constant at 10''^ erg, for different values of n. In the lower panel, the fraction of destroyed dust mass is plotted (for different n) as a function 
of explosion energy E, keeping the ejecta mass constant at 10^'' g. 



